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Abstract 

Vortex breaking has been traditionally studied for nonuniform critical current densi- 
ties, although it may also appear due to nonuniform pinning force distributions. In this 
article we study the case of a high-pinning/low-pinning/high-pinning layered structure. 
We have developed an elastic model for describing the deformation of a vortex in these 
systems in the presence of a uniform transport current density J for any arbitrary ori- 
entation of the transport current and the magnetic field. If J is above a certain critical 
value, J c , the vortex breaks and a finite effective resistance appears. Our model can be 
applied to some experimental configurations where vortex breaking naturally exists. This 
is the case for YBa2Cu30y_,5 (YBCO) low angle grain boundaries and films on vicinal 
substrates, where the breaking is experienced by Abrikosov-Josephson vortices (AJV) 
and Josephson string vortices (SV), respectively. With our model, we have experimen- 
tally extracted some intrinsic parameters of the AJV and SV, such as the line tension e; 
and compared it to existing predictions based on the vortex structure. 
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1 Introduction 

The study of vortex physics has been an important topic in superconductivity research since 
the prediction of the existance of superconducting vortices by Abrikosov [Tj [2j [31 0]. Su- 
perconducting vortices have been the subject of many theoretical studies including those 
addressing the vortex line energy [5j [6] , single vortex interactions EJ QH] and vortex 
lattice interactions [TT1 H2l [T3 | Hi] . 

One way to measure the properties of a vortex, such as the vortex line tension ei is by 
means of vortex breaking experiments. These kind of experiments were initially performed 
by using a nonuniform transport current in order to create a nonuniform driving force in the 
vortex length thus bending the vortex and, for large enough current density difference, break 
(or cut) the vortex [15| . I16 |. ITT ^. I19 |. I18| . Another way in which vortex breaking may be observed 
is where vortices are subject to an inhomogenous pinning force in the presence of a uniform 
transport current. In this article, we study the microscopic bending and breaking process for 
the situation where the vortex crosses a low-pinning region (LPR) of thickness d compressed 
by two high-pinning ones (HPR) with a uniform transport current J and a magnetic field 
B with an arbitrary orientation, Fig. [TJ In order to perform this study, an elastic model is 
developed assuming isotropic vortices with a constant line tension ei towards deformations. 
Using this model, we propose a vortex breaking process in the depressed pinning region. 
Afterwards, a J c (9,ip) dependence is obtained, where 9 and (p are the angles defining the 
orientation of the applied field (see Fig. [1]), as a function of the elastic force per unit length 
and the pinning line forces in either the LPR and the HPR. 




Figure 1: (a) Sketch of an undcformed vortex line in a depressed pinning layer, consisting of a low 
pinning region (LPR) bounded by high pinning ones (HPR). The applied field B Q makes an angle 
with both the LPR-HPR interface and J. The vortex length in the LPR, D, is D — d/\six\6 svmp\. 
(b) Definition of the angular coordinates (p and 8. 
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The geometry described in this article can be directly applied in order to model experimen- 
tal situations like low-angle grain boundaries (LAGB) in YBa2Cu30y_5 (YBCO) bicrystals or 
YBCO films on a miscut substrate (commonly named vicinal films), as discussed in Sec. [U It 
is worth noting that, according to Gurevich et al [20] the vortices in high-temperature super- 
conductor LAGB for crystal misorientation angles up to 23° are not conventional Abrikosov 
vortices (AV) but Abrikosov-Josephson vortices (AJV)0. Moreover, the vortices that it is 
possible to break in YBCO vicinal films are string vortices (SV) lying between the ab planes, 
which are of Josephson nature. In this paper, the vortex line tension e\ in these systems 
is found by fitting the theoretically obtained J c (8,(p) dependence to measured data. Then, 
thanks to the model developed in this work, it is possible to extract the line tension of exotic 
superconductors of Josephson and Abrikosov-Josephson nature, which physical properties are 
still not very well known. 

The article is structured as follows. First, some general features of our problem, including 
the basis of the elastic model, are discussed in Sec. EJ In Sec. [3] we deduce the analytical 
limits of the vortex deformation and critical current density for the studied HPR/LPR/HPR 
structure at two specific magnetic field orientations: B perpendicular or parallel to J, with J 
crossing the HPR/LPR boundary. In the following section, Sec. 01 we introduce a numerical 
procedure for solving the problem for any field and current orientation and we present general 
normalized results. Next, Sec. [5j the model is applied to the experimental situations of LAGB 
in YBCO bicrystals and vicinal films in order to obtain the AJV and SV vortex line tension, 
respectively. In that section, it is also discussed the flux breaking process in a vicinal film. 
Finally, we present our conclusions in Sec. 



2 General considerations 

In this section we discuss some general features of our system, such as the basis of the elastic 
model and the properties of the vortex line energy. 



2.1 Elastic model 

We simulate the behaviour of a vortex line as an string with constant line energy e/. Thus, 
its total energy is 



E = ei dl = ei dr\d T r\, (1) 
Jc Jc 

where C is the vortex line path, / is the line length and r is a parametrization of C. If the 
vortex follows a straight line, the elastic force at the vortex ends is 

F e = ±e ; d£/d/ = ±eie l: (2) 

where e/ is the unit vector locally parallel to the vortex line: 

e,M = (3) 

From Eq. (2) it is clear that for an isotropic string the line tension equals to the line 
energy. The local elastic force per unit length f e (r) for a deformed vortex can be deduced by 
approximating the vortex line into a set of straight segments and using that the net elastic 
force at a junction between segments is 

SF e (r) = ei[ei(r) - e t (T - St)], (4) 



Actually, the structure of a vortex in a LAGB for anisotropic superconductors is more complex than the 
AJV described in 1201. The vortex nature in our studied situation is discussed in Sec. 15.11 
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where 5t is a variation in r and SF e (r) is the elastic force in the junction between the vortex 
segments corresponding to r and r — 5t. After making the limit of 5t — > 0, the force per unit 
length can be found as 

f e (r) = j^d T r x (d 2 T r x d T r). (5) 

Apart from the elastic force, the vortex is submitted to a driving force, F^, due to a 
macroscopic current density J. The line force density of this force, f^, is 

f d = $ J x ej, (6) 

where $o is the flux quantum, that in SI is <&o = ^r^e. 

For superconductors with pinning we should also take into account the pinning force 
per unit length f p , which has a certain maximum magnitude f pm - This force opposes any 
local movement of the vortex line, so that it will follow the direction of the elastic and the 
driving forces. This description of the pinning behaviour implicitly assumes that f p is due to 
many weak pinning centres uniformly distributed in the vortex length, rather than few strong 
pinning sites. The latter scenario would result in a more complicated bending behaviour. 

The equilibrium vortex configuration will be that one which results in a null net force per 
unit length at any point. Thus, from Eqs. (|5|) and ([6]) we find the differential equation of the 
vortex line 

T^d T v x (^r x d T r) + $ J x + f p = 0. (7) 



2.2 On the line tension 

We next discuss the dependence of the line tension e\ on the superconductor internal param- 
eters for both an isotropic and anisotropic Abrikosov vortex (AV). 

For an isolated isotropic AV, the line tension corresponds to the energy per unit length, 
which has been investigated by several authors EJ E] . Following the numerical calculations 
of Hu [6j , ei in the approximation of high Ginzburg-Landau parameter k is 

e, = e,, (In as + 0.4968), (8) 

where e^o = (^ > o/^) 2 /(47r//o) and A is the magnetic penetration depth. A more general, 
although maybe less accurate, analytical expression for arbitrary k was obtained by Clem [5] 
using a variational approach. That expression, however, is still proportional to ej o- 

A vortex in a flux line lattice has a different line tension. In general, q depends on 
the vortex separation, oq, and the characteristic deformation wavelength, [3J 0]. The 
qualitative behaviour is that t\ decreases with decreasing both ao and Id- However, for 
magnetic fields such that oo < A and ao 3> £, the line tension is roughly constant (where £ is 
the superconductor coherence length). In a simplified way, this can be explained as follows: 
if ao <C A the magnetic field and the supercurrent density is roughly uniform between the 
vortices except close to the vortex core. Then, the energy variation per unit length due to a 
deformation is just the core energy per unit length and, consequently, it does not depend on 
the magnetic field. Following [10], e/ for the London limit is ei = e^o dx(e~ x ^ K )/Vx 2 — 1 + 
e«,c = Q,o-^o(l/^) + Q,c! where e; jC is the core energy per unit length and Kq is the modified 
Bessel function of order 0. Comparing with Eq. ([H]), we obtain that 

e,, c = e/,o [In k + 0.4968 - K (1/k)] . (9) 

For the limit of k > 1, Eq. ([9]) become^] e\ c = 0.38096/^, so that we obtain a line tension of 
the order of e^o. 

2 The approximation of k S> 1 produces an error smaller than a 1% for k > 20. 
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In Sec. [5]we study the line tension in SVs in YBCO, for which q is essentially the same as 
an anisotropic Abrikosov vortex [31 13] • Following Ref . [3] , the line tension for an anisotropic 
vortex parallel to the ab planes with k 3> 1 and El Id *C \ a b, where X a b is the penetration 
depth for the ab planes, is approximately 

*t - r-^V (io) 



A 2 



r 47r^ A, 



(11) 

ab 



where ej~ and ej are for deformations in the c and ab directions, respectively, and T is the 
material anisotropy factor T = X c /X a b > 1 (for YBCO, F 5), where A c is the penetration 
depth in the c axis. 



3 Analytical limits 

In this section we analytically describe the limits of (p — ► with 9 = tt/2 (B perpendicular 
to J) and <p — > tt/2 with 9 = tt/2 (B parallel to J). 



3.1 Applied field perpendicular to the current flow 

In the following, we suppose that B is almost parallel to the HPR-LPR interface, that is very 
low ip and 9 = tt/2, Fig. [U For this case, the applied field is approximately perpendicular to 
the current flow, so that the vortex bows in the z direction only. This is consistent with the 
numerical calculations in Sec. 14.31 below. 

Let define the axis x' as that in the direction of the applied field and y' as that in the 
direction of e z x e x i, Fig. [TJ Taking the coordinate x' as the vortex line parameter, its 
positions are 

r(x') = x'e x , + z{x')e z , (12) 

so that the z(x') dependence uniquely describes the vortex line. Using this parametrization, 
it can be seen from Eq. ([5]) that the elastic force per unit length is in the e y / x e/ direction 
and, thus, it can only compensate the driving force in this direction. Therefore, the relevant 
component of the driving force is ^ofJ • e y '] e y' x e z = &oJ cos ¥ e y' x e i- The differential 
equation for the vortex line of Eq. (JT]) becomes 



J<&o cos (p — €i 



dz 
cb 7 



"3/2 ,2 



^ + / P = 0. (13) 



3.1.1 Infinite pinning force in the high pinning region 

In order to simplify the system, we assume that the maximum pinning force per unit length 
in the high pinning region is much higher than any other line forces in the problem, i. e. the 
pinning force in the low pinning region, the driving force and the elastic force. With this 
approximation, the vortex in the HPR follows a straight line in the B direction for any J. 
Thus, we only have to obtain the vortex line configuration in the LPR. First, we study the 
case of null pinning force in the LPR and afterwards we include the effect of a finite / Pi m- 

3 For our experimentally studied situations of Sec. the condition Id \ a t, is satisfied for ip > 5° and 
8 > 0.5° for LAGB and vicinal films, respectively. 
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Assuming null pinning force, the differential equation (1131) can be solved by direct inte- 

dz 
Ax' 



gration. For doing this, Eq. (| 13j) must be first solved for 4^t(x') and then for z(x') with a 



result 



z(x') = h- ^R 2 -(x' + k 2 ) 2 , (14) 

where R is defined as R = e//(| cos<^| J^o) an d k\ and k 2 are integration constants. These 
constants can be found taking into account that the vortex line in the HPR is fixed, obtaining 
that 

z{x') = ^R 2 - (D/2) 2 - VR 2 - x' 2 , (15) 

where D = d/\ sin#sin<^| is the length of a straight vortex in the LPR (Fig. [JJ). From Eq. 
(]15p we see that the vortex line bents forming a circular arch with radius R. The elastic and 
driving forces per unit length are calculated from Eqs. ©, © and (fT5"|) . yielding the result 

f e = --^ e y' x e l (16) 
i d = cos (p$ Je y/ x e;. (17) 

^From these equations we see that f e and have uniform magnitude and are perpendicular 
to the vortex line. The elastic force per unit length is inversely proportional to the curvature 
radius R. Therefore, the maximum |f e | that the vortex is able to produce, f e m, is 



f e>m = 2ei/D (18) 



that leads to a critical J, J c , as 



J c = ^r^ = 1% tan \<p\, (19) 

Where at the second equality we used that D = d/ sin \ip\ for 8 = n/2 (Fig. [T]). For J lower 
than J c the vortex line is stable, while for J > J c it cannot balance the driving force and 
the vortex line breaks, Fig. [2j For J > J c , Fig. 0c), there is a straight horizontal vortex 
segment in the border between the HPR and LPR because in this place the driving force 
pushes the vortex towards the HPR, where the available pinning force is very large and, thus, 
fully compensates the driving force. 

The breaking process is completed if the horizontal vortex segment cross-joins and re- 
combines [TJ EJ [10] with the fixed part of the neighbouring vortex. This will happen when 
the driving force overcomes the repulsion between these two segments. Indeed, not only is 
the repulsion between vortices at right angles strongly suppressed but also the interaction 
force becomes attractive for separations close to £ for any orientation [7J [9]. The expected 
situation after cross-joining is that the vortex segment in the LPR joins with the undeformed 
neighbouring one and there appear two horizontal antiparallel vortex lines. These horizon- 
tal vortices will experience a strong attraction which can overcome the pinning and driving 
forces, so that the vortices move towards each other and anihilate. In the case that the 
vortex attraction is not strong enough, the elongation and recombination process repeats, 
creating horizontal vortices with multiple vorticity. These vortices have a larger core and, 
thus, present a lower pinning line force and their mutual attraction is larger. At some mo- 
ment these horizontal vortices will detach from the HPR-LPR border and anihilate to each 
other. Afterwards, the breaking and cross-joining process is repeated. 

For large enough magnetic fields, i. e. low vortex separations, there can be vortex cross- 
joining for J below that of maximum vortex deformation when the distance between the 
vortex arc in the LPR and the straigh part of the next vortex is of the order of £. 
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Figure 2: Vortex bowing for the situation of Fig. [T]with = 7r/2, ip <C it/2 and infinite maximum 
pinning force per unit length in the HPR. The graphs correspond to J < J c (a), J — J c (b) and 
J > J c (c) from left to right. 



Let define the vortex breaking force Fbreak as the required driving force magnitude nec- 
essary to break the vortex in absence of pinning in the LPR. Using Eqs. ([3]) and (|6]) we find 
that the driving force is 



f f D/2 dx' 

JLPR J-D/2 



r 



/LPR J-D/2 \d x >r\ 

,D/2 

= cos ip&oJ e y x / dx d x /r = — cos tp&o J De z . (20) 

J-D/2 

The breaking force is obtained by taking J = J c and using Eqs. (fT9j) and (f20|) . with a result 

Fbreak = 2e/. (21) 

We next consider the presence of a pinning force in the LPR as follows. If the driving 
force per unit length is lower than / Pim , the pinning force per unit length counteracts the 
driving force and the vortex stays as a straight line. Otherwise, the vortex bows and the 
elastic force per unit length compensates the remaining force. For this situation, the vortex 
line follows Eq. (|13j) with f p = —fp. m . Taking this into account, the vortex configuration still 
follows Eqs. (fTij) and (fT5|) but with a curvature radius R = e;/(| cos <p\ J<I>o — fp,m)- Thus, 
Eqs. (fT6l) - (fT8|) are still valid but now J c is 

Je = jjfL tanM+ (22 ) 
$o|cos</?| 

This equation is equivalent to that deduced by Durrell et al in [211 [22] . if it is taken into 
account that Fb rea k = 2e; . 

The vortex line of Fig. [2] is nonderivable at the border of the LPR and HPR. Then, on 
these points there acts a net elastic force with magnitude |F e | = ei\e x ' — e/(x' = d/2)\ [Eq. 
(J3j)] which can only be balanced by a finite pinning force. Such a pinning force is possible 
under the assumption of infinite maximum pinning force per unit length in the HPR, fp m . 
If f* m is limited to a finite value, the vortex line must be derivative at every point. Thus, 
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the vortex can bow in the HPR even if the driving force per unit length is lower than f pm . 
However, the main features of the vortex breaking remain the same. 

We contemplate a finite pinning force per unit length in the HPR in the following section. 

3.1.2 Finite pinning force in the high pinning region 

We next relax the assumption of infinite maximum pinning force per unit length in the HPR, 

J p,m' 

The vortex configuration z(x') still follows the differential equation (I13p but now f p de- 
pends on the position. For f* m > QqJ > f P;m , the vortex line in the LPR bows in the 
direction of the driving force and, then, f p = —f P) m there. The pinning force per unit length 
in the HPR can, in principle, take any value between —f£ m and /p >m . If the transport current 
density is reached by increasing its value from to J in a quasistatic speed, the vortex line 
progressively deforms from a straight line to its final stable configuration. In this process 
the pinning force density opposes with its maximum value to the driving force. Thus, in the 
HPR fp = f pm . 

Taking this into account, we solve Eq. (|13|) using that the vortex line must be continuous 
and smooth at the LPR and HPR boundaries, with the result 



z(x') 



for x' < —x' 

R* 2 - (x' + x' ) 2 - R* for -x' < x' < -D/2 

^R 2 - x' 2 + {2x' /D)^R 2 - D' 2 /4 - R* for - D/2 < x' < D/2 (23) 

R* 2 - (x> - x' ) 2 - R* for D/2 < x' < x' 

for x'q < x' 



with 



R = e { /(|coe ¥ >|J* -/p,m) (24) 
i?* = e / /(/* OT -|cos ¥3 |J$ ) (25) 
x' = (D/2)(l+R*/R). (26) 

The quantity x' is the minimum \x'\ for which the vortex is not deformed and R and R* are 
the curvature radii in the LPR and the HPR, respectively. 

The vortex breaking process for this situation is shown in Figs. [3]and[31 where we took 
fp,m/fp,m = 9 and 2ei/(df Pjm ) = 18.6, corresponding to the experimental situation in Fig. 
El From Figs. [3] and 0] we can see that the breaking process is essentially the same as 
for infinite f* m (Fig. However, for J below that of maximum elastic force, the vortex 
bending could cause a local separation between neighboring vortices smaller than £. This 
would reduce J c compared with the predicted one from elasticity considerations. For our 
experimental situation in Sec. 15. Q the minimum vortex distance at J = J c is around 2£ so 
that the expected cross-joining is still at the breaking critical current density. 

For local vortex separations larger than £, such as those in Figs. [3] and H] the limiting 
J c mechanism is basically the same as for infinite /* (Fig. [2]). Thus, Eq. (fT8|) for the 
maximum elastic line force is still valid. However, in order to observe breaking for finite f p ^ m , 
it is not enough that fj overcomes f p>m but also 

fp,m > fp,m + fe,m = f P ,m + 2e z | sinip\/d. (27) 

This can be seen as follows. The driving force per unit length is uniform in the vortex which, 
from the equilibrium condition in the LPR it is fd = f p ^ m + f e , m - Then, if f* m < fd = 
fp,m + fe,m, there is depinning in the HPR and there is no vortex breaking. Moreover, from 
equation (|27"|) it can be seen that 2ei/d > f* m — f Ptm there will be breaking or depinning in 
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Figure 3: Vortex bowing for the situation of Fig. [T]with 8 = tt/2, <p = 10° and a finite maximum 
pinning force in both the LPR and HPR f Pim and /* m , respectively. The graphs are for f* m / f P , m = 9, 
2ei/(df p>m ) — 18.6, d — lOnm and B = 6T, corresponding to the experimental situation in Fig. [H] 
Half the number of vortices have been plot for a better visualization. 



the HPR, depending on the value of ip. Thus, from Eqs. ([2"2"]) and ([27]) . the J c (p) dependence 
will be 

f Jc = #77 tan Id + a { p - m I for <p < ip d 
1 for^>^ 

with 

<p d = arcsm ^ . (29) 

Another difference with the assumption of infinite /p TO is that now the vortex significantly 
bends in the HPR, with a maximum deformation at J = J c . The curvature radius in the 
HPR at J = J c , R*, can be deduced from Eqs. (I25p and (22), with a result 

R*/d = . (30) 

fp, m -f P ,m- -^/l sm^| 

The value of R*/d is finite for the limit y> = and R*/d approaches to infinite when \ip\ = p d . 
The behaviour of R*/d as a function of p is plotted in Fig. [5] for 2ei/(df Pim ) = 19 and 
fp,m/fp,m = 5,10,15,20,25 and 30 in the arrow direction. The angle ip^, calculated from 
Eq. J2SI), is 12.2,28.3,47.5 and 90° for f* m /f p , m = 5,10,15 and 20, respectively. In Fig. E 
it can be seen that R* diverges when it approaches to <p> d but it has no discontinuities for 
2ez/(rffp,m.) < fp,ml fp,m ~ 1- 

3.2 Applied field parallel to the current flow 

We next describe the situation that the applied field is parallel to the current flow, usually 
known as the force-free configuration. For simplicity, we only study the case of infinite pinning 
force per unit length in the HPR and no pinning in the LPR. 
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Figure 4: The same as Fig. [3] but for ip = 20°. 



This geometry corresponds to(£> = # = 7r/2in Fig. [H so that B = Be y . Taking y as the 
vortex line parameter, the differential equation (|7|) turns into 



(d%xe x + d%ze z ) $ J 



H (9 y 2;e x - fl^xe*) = 0. 



(31) 



For the force-free configuration, it has been found that the vortices are susceptible to 
become spirals |23l I24| . Moreover, our numerical results below (Fig. m Sec. H]) show that 
for a B almost parallel to J, the vortex line configuration is an spiral with twist pitch —d. 
The spiral must cross the y axis at y = ±d/2 for continuity with the vortex line in the HPR. 
A spiral following these features is 



r( U) = a [ 1 + cos I e, 



2iry 

ye v + asm — — e z 

a 



(32) 



where a is the radius of the spiral. This spiral follows Eq. (|31|) and, thus, it describes a static 
vortex configuration. The spiral of Eq. (|32p has the axis in the y direction but displaced a 
vector ae x form it. Actually, any rotation of the spiral from Eq. (|32l) around the y direction 
is also a static vortex configuration. The elastic and driving forces per unit length for these 
spiral vortex lines are 



i + (tt) 



2aiT 
d 



--e t 



with 



cos 



( 2wy 
V d 



e x + sin 



2vry 



(33) 
(34) 

(35) 
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10 20 30 40 50 60 70 80 90 

<P(°) 

Figure 5: Vortex line radius in the HPR for J = J c normalized to d, R*/d, as a function of ip 
(Fig. E) for 2ei/(df p>m ) = 19 and f* im /f p , m = 5, 10, 15, 20, 25, 30 in the arrow direction. The line for 
fp, m /f P ,m = 20 diverges for (p = 90°'. 



The radius at which there is equilibrium a eq is deduced from Eqs. ([55]) and and consid- 
ering f e + frf = 0, obtaining 



J^oJ V2vr 



(36) 



By examination of Eqs. (|33|) and (|34|) we see that the vortex at a = a eq is in unstable 
equilibrium. Indeed, a radius slightly smaller than o eq presents an elastic force higher than 
the driving force and the spiral shrinks until it becomes an straight line. For a > o eq the 
driving force is higher than the elastic one and the vortex expands indefinitely. However, any 
equilibrium is not possible when J is high enough to produce an a eq from Eq. f)36[) with null 
or imaginary values. The minimum J value at which this occurs, J c , is 



Comparing J c for the parallel and the perpendicular cases, Eqs. (|T9|) and (j37j) . we see 
that for the same vortex length in the LPR J c for the force-free situation is only tt times 
higher than that one for a perpendicular field. 



4 Numerical approach 

In this section a numerical model is introduced in order to calculate the vortex deformation 
in a depressed pinning layer, Fig. [H at any orientation of the macroscopic current density 
J and the applied magnetic field B. As for the analytical approximation, we consider an 
elastic model for the vortex line. In order to simplify the analysis, we assume the approx- 
imation of infinite maximum pinning force per unit length in the HPR, f* m — > oo, for the 
numerical calculations. As shown above, Sec. 13.1.21 taking into account a finite f* m does 
not substantially change the J c results. 
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4.1 Vortex line calculation 



We next outline the numerical method used for calculating the vortex line configuration. In 
order to generalize the procedure for several applications, we consider a vortex segment of 
length d strongly pinned at its ends, beyond them the vortex line is straight and follows the 
applied field direction. 

The vortex configuration is calculated by physical evolution in a similar way as in [TU] , as 
follows. 

We take the initial state of a straight vortex in the applied field direction and we divide 
the pinning-free segment into n identical line elements. In the following discussion we term 
the boundary between elements as a node. Then, we calculate the forces per unit length f e , 
id and i p corresponding to each node and change the node position proportionally to the 
resulting line force. That is, the position of the node labelled as i at the time instant k, r^, 
is calculated from the previous time instant as 

*i,k = + A fc [f e (ri )fc _i) + fd(r i)fc _i) + f p (r i)fe _i)] (38) 

with 

A k = (t k - t fc _i)/i7, (39) 

where r\ is a proportionality factor and t k is the time at instant k. The time evolution of Eqs. 
(|38p -(|39 p has the physical basis of an asymptotic viscous vortex movement with a viscosity 
per unit length rj, for which r/r = i e + id + i p |1Q|. [3] . In the numerical method, we implemented 
a self-adaptative procedure for choosing the optimum at each time instant for a quick 
convergence to the stationary situation. 

The elastic and driving line forces, f e and id are calculated using Eqs. (JH) and (J6j) , 
respectively, and i p is evaluated taking into account that it opposes to the action of the 
driving force and has a maximum magnitude f p , m - Then, the effect of i p when the driving 
line force is smaller than f p>m is to avoid any vortex deformation for |f^| < / Pjjri and for higher 
id, i p reduces the effect of the driving force by / p , m . 

4.2 Critical current calculation 

For a vortex breaking phenomenon, the critical current J c is the maximum J for which there 
exists a stable vortex line configuration. A current above J c creates a driving force per unit 
length that cannot be balanced by neither the elastic nor the pinning ones and the vortex 
moves or deforms indefinitely. 

According to this, the critical current is calculated as follows. We take a J interval with 
one boundary below J c , J m \ n , and the other one above, J max - Then, we choose a certain 
current between them, J', and calculate the vortex distribution. If the vortex line for J = J' 
converges to a stable solution, we change J m \ n into J' and J max into J', otherwise. Thus, 
the interval containing J c is narrowed. In order to ensure a solution of the procedure, we 
choose a very wide initial interval, with J m i n = 10 3 A/m 2 and J max = 10 15 A/m 2 . We found 
that the number of vortex line calculations is minimized by taking J' as the middle point in 
logarithmic scale, J' = V Jmin^max- The number of vortex line evaluations for obtaining a 
1% wide interval is around 12. 

Numerical calculations showed that the J c result do not significantly change with the 
number of vortex elements n. The computing time for n = 20 in an standard table computer 
is of few seconds (around 15s), allowing a systematic study of the J c dependence with the 
applied magnetic field orientation. 
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4.3 Results 



4.3.1 Vortex bending 

Although the numerical method described above is valid for any relative orientation between 
J and B, here we restrict our study for B in the xy plane of Fig. [U that is, 9 = tt/2. This 
corresponds to the experimental situation of a thin film YBCO bicrystal with the applied 
field parallel to the surface and, consequently, the ab planes. 

In Fig. [6] we present the vortex bending at J = J c for several p orientations. In the plots, 
the x' axis is in the direction of the applied magnetic field and the y' axis is perpendicular 
to both x' and z, Fig. [TJ 

The situation in Fig. [6^a) is for J almost perpendicular to the field, with <p = 0.027r(3.6°). 
It can be seen that the vortex line bends almost exclusively in the z direction, following a half 
circle with a diameter equal to the vortex length in the LPR, D, as predicted analytically in 
Sec. 13.11 When J makes an intermediate angle with B, as in Fig. E^b), the vortex also bends 
antisymmetrically in the y' direction and the maximum deformation in the z axis is smaller 
than D. Another interesting situation is when J is almost parallel to B, as shown in Fig. 
d^c) for which cp = 0.497r(88.2 o ). For this limit, the vortex approaches to an helix with twist 
pitch —D. This result justifies the helical vortex assumption for the analytical deduction in 
Sec. EE2J 



4.3.2 Critical current density 

In the following we study first J c with null pinning in the LPR and then we consider the 
effect of a finite f Pt m- Below, we refer to J c for f prn = as the breaking current density, 

Jbreak • 

In order to study the general breaking behaviour in a GB, in Fig. [7J (solid line plus 
symbols) we present Jbreak normalized to the breaking current density at the force-free sit- 
uation Jbreak,ff = 2"7re//((f<3?o) [Eq. (f3"T|) ] as a function of the angle ip (see sketch in Fig. [T]). 
Such representation is universal for any d and e\. This can be seen as follows. Jbreak is the 
maximum J at which the elastic line force compensates the driving one. Since the elastic 
line force is proportional to q, the same must be for Jbreak- Besides, the whole vortex line 
configuration at J = J c scales with the LPR thickness d. Since f e is proportional to the 
vortex line curvature and twisting, f e and, thus, Jbreak must be inversely proportional to d. 
Consequently, Jbreak is proportional to ei/d and to Jbreak.ff = 27re// (cZ^>o)- 

The maximum Jbreak (Fig. [7j) is for <p = tt/2, corresponding to J parallel to B. As can be 
seen in the figure, the Jbreak value at the maximum is Jbreak.ff; i n accordance to the analytical 
deduction in Sec. 13.11 

The numerically calculated Jbreak /Jbreak,ff of Fig. [7] can be well fitted by the expression 

1 0.22 

Jbreak/ J c ,// = - tan |0.78v?| + - — sin \tp\, (40) 

TT TT 

corresponding to the dash line in Fig. [JJ The difference between the analytical fit and the 
numerical data decreases with decreasing ip, being of a 5% for p = tt/2 and less than 1.3% 
for \p\ < 0.4vr(72°). 

We next study the effect of including a pinning force in the LPR. In Fig. [7J we plot 
Jc/Jcjf for a nonzero f P:m [fp,m/(<&oJcjf) = 0.02,0.1], showing a roughly uniform increase 
of J c by fp-m/^Q in most of the J c (p) curve except close to (p = tt/2, where J c diverges. 
This behaviour can be explained as follows. If / Pjm is large enough to counteract the driving 
force, the vortex stays as a straight line. Otherwise, the vortex bends and the elastic force 
adds to the pinning one. However, the deformation of the vortex does not always lead to 
a compensation of the driving force. Indeed, for J almost parallel to B and high enough 
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Figure 6: Vortex deformation for the situation of Fig. [T]at J = J c for — ir/2 and several <p angles 
(red thick lines), where all dimensions have been normalized to the GB half width d/2. Blue thin 
lines are the projections in the x'z, x'y' and y'z planes, the black arrow shows the current direction 
and the grey lines show the LPR-HPR boundaries in the x'y' plane. In (c), note the different scales 
of the axis. i A 




9 

Figure 7: Calculated critical current density, J c , for 8 = ir/2 as a function of <£> for the situation of 
Fig. [1] The curves are universal, for any d and e; when, J c is normalized to the critical current density 
in the force-free configuration Jbreak,ff, expressed by Eq. ([57)) . Lines plus symbols are for numerically 
calculated data for / p , m /(^o Jbreak,ff) =0,0.02,0.1 from bottom to top and the dash red line is for the 
analytical fit of Eq. 
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J, the vortex deforms helically and experiences a higher driving force than when it was 
straight. Then, for this situation the critical current density is determined by depinning, 
so that J c = f p>m /($a\ cos <p\) and J c diverges for tp = tt/2. For the orientations that flux 
breaking exists, the increase in the numerically calculated J c due to pinning in the LPR fits 
well to /p, m / (^ovl cos ip\). Thus, taking into account the above features and Eq. (jlOl) for 
-Tbreak) Jc can be well described by 

J c = max(J C) i, J C)2 ) (41) 

with 

= ifei < 42 > 

J c,2 = (tan 10.78^1 +0.22 sin + —Js&—. (43) 

d®o 3>ov|cos(p| 

Comparing Eq. (|43|) and equation (|22|) for low magnetic field orientations p, it can be seen 
that for the low-</? limit they are coincident and they differ less than 2% and 8% for ip < 15° 
and (p < 30°, respectively. 

If we take into account that the pinning line force in the HPR has a finite value, the critical 
current density cannot be larger than that limited by depinning in the HPR (Sec. 13. 1 .2j) . 
Thus, J c will have the lower value between those given by Eqs. (|4*T|) - (j4*3|) and f^ m / (<&o| cos p\). 



5 Application to experimental situations 

In this section we apply the theoretical model developed in Sees. [3] and H] to the experimental 
situations of LAGB and vicinal films in YBCO. The deduced J c (0, p) is fitted to the measured 
data in order to obtain the vortex line tension in these systems. 



5.1 Low-angle grain boundaries 

The study of low-angle grain boundaries is of great practical importance because they are 
present in YBCO coated conductors [25], which are the most promising high-temperature 
superconducting materials for applications in electrical devices [26J. Indeed, grain boundaries 
are limiting the critical current density in coated conductors. 

The geometry described in this article can be directly applied in order to model low- 
angle grain boundaries in YBCO bicrystals. In these grain boundaries there are Abrikosov- 
Josephson vortices (A JV) [20J , which magnetic extent and normal core is of larger size than 
the usual Abrikosov votices and, thus, the AJV in the grain boundary presents lower pinning 
than the AV in the grains [201 [27] . Assuming a sharp transition between the grain and the 
grain boundary, this situation corresponds to the HPR/LPR/HPR layered structure in Fig. [TJ 
Actually, for the Abrikosov-Josephson vortices in LAGB introduced in Ref. [20] it is assumed 
that the vortices in the crystals are conventional AV and parallel to the boundary plane. The 
AJV theory in Ref. [20] is, then, directly applicable to high-temperature superconductors for 
magnetic inductions and LAGB parallel to the c crystallographic axis. For B in the ab plane, 
the vortices in the crystals are not AV but Josephson string vortices (SV), which complicates 
the description near a LAGB. However, it is expected that SV in a LAGB will present a 
larger phase core and magnetic extent in the LAGB direction. This significantly changes the 
SV properties, as follows. 

The supercurrents around SVs flow in two characteristic penetration depths, X a b in the c 
axis direction and A c = TA a b in the ab plane 011]. In SVs, the phase core have a size of dn 
in the c axis and Td[\ in the ab planes, where d[\ is the separation between ab layers. SV are, 
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then, highly anisotropic. A LAGB perpendicular to the SV and parallel to the c axis will 
enlarge both the penetration depth and the phase core in the c direction, reducing the SV 
anisotropy. Furthermore, SVs present a strong intrinsic pinning towards driving forces in the 
c direction owing to the reduced order parameter between ab planes. Therefore, increasing 
the SV phase core in the c direction due to a LAGB will significantly reduce the intrinsic 
pinning. If the SVs are making a certain angle ip with the LAGB, the situation is even more 
complex but for not very large <p it is expected a significant intrinsic pinning reduction with 
a weak ip dependence. 

5.1.1 Comparison with experiments 




Figure 8: Measured J c for a grain boundary in an YBCO bicrystal with a 4.9° misoricntation angle [^Tj 
(hollow blue circles) together with the theoretical fit for breaking in the LAGB (solid line), calculated 
from Eqs. (@T])-(jin]), and for depinning in the crystals (dash line), from J c = /p, m /($o| cos^|). The 
fit parameters are e ; = 2.89 x KT 13 N, / ftTO = 3.1 x 10~ 6 N/m and f* m = 2.8 x HT 5 N/m. 



The J c numerical results from Fig. \7\ (Sec. 14. 3p . or the equivalent formulae (l4*Tj) - (j4*3l) . 

can be used for fitting experimental data for grain boundaries and extract the vortex line 
tension In Fig. [5] we fit the experimental J c (f) dependence for a low angle GB in an 
YBCO bicrystal |21j with our numerical curve, obtaining e, = F break /2 = 2.89 x 1(T 13 N, 
f P ,m = 3.1 x 10~ 6 N/m and /* = 2.8 x lCT 5 N/m for a temperature T = 65K and an applied 
magnetic field B = 6T. From the fitting parameters and Eq. (|29H a maximum flux breaking 
angle of ip^ = 25.5° is obtained, corresponding to the measured value in Fig. [HJ Moreover, 
the experimental results from [21] are field independent in the measured range (1 to 5T). 
This can be explained by the low vortex separation compared to \ a b, which causes that the 
line tension is mainly the vortex core energy, as discussed in Sec. 12.21 

It is interesting to compare the fitted value of ei with the expected one for an anisotropic 
vortex. Using a X a b temperature dependence of 

Kb = \abfi{l-T/T c )- 1 ' 2 (44) 

with A a{ , i0 = 1.50 x lCT 7 m and T c = 91K, Eq. (JTOD yields e t 1.72 x KT n N, which is two 
orders of magnitude larger than the value obtained from the experimental data. However, a 
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better agreement is found when assuming an isotropic vortex with a magnetic size A c = rA a 6, 
obtaining e/ ~ e^o = 1-38 x 10 _13 N. This is consistent with a significant vortex anisotropy 
loss due to a vortex size increase in the c direction, caused by the LAGB (Sec. I5.ip . However, 
the origin of this discrepancy might be the oversimplified analysis in the deduction of Eqs. 
(I10p - (llip for the line tension. A more accurate prediction would require a further complicated 
discussion on the vortex lattice elasticity [T3J H] ■ 

5.2 YBCO films with miscut substrate 

Another system which spontaneously presents flux breaking is an YBCO film on a miscut 
substrate [281 [22], commonly named a vicinal film. It is well known that for applied fields 
with an orientation angle with the ab planes below around 11°, the vortex line in YBCO have 
a stair-like shape, composed by pancake vortices (PV) in the ab planes and string vortices 
(SV) between these planes [3J. While pancake vortices are strongly pinned, this is not the 
case for the string ones. Although a large pinning force is required to depin SVs from the ab 
plane direction, they are almost free of pinning for forces parallel to the ab planes. The result 
is that the SV is a weakly pinned vortex segment confined between two ab planes with its ends 
strongly pinned by pancake vortices. Therefore, under the presence of a transport current 
with a component in the c direction, the SVs will bend in the ab plane and, eventually, break 
for a transport current density above a certain threshold^!, J c . 



n 




>j 



a 



Figure 9: Sketch of a vicinal film and description of the angular coordinates 9 and ip. 



5.2.1 Breaking description for vicinal films 

Due to intrinsic pinning, the driving force will deform the SV mainly in the ab plane following 
a circular arch, as shown in Fig. [TOl a.b). If we restrict our analysis to the vortex deformation 
in the ab planes only, a SV would expand until it finds another string vortex in the same plane. 
However, flux crossing will occur much sooner if we consider the interaction between SV in 
different planes. For J = J c , Fig. fTOT b) . the SVs between different planes are antiparallel 
at the junction with the common pancake vortex, where they experience a strong attractive 
interaction. This attraction compensates, at least partially, the bowing effect of the driving 
force and antiparallel straight SV segments appear [TOlch When the length of these segments 
is large enough, the SV attraction will overcome the intrinsic pinning and they will cross-join, 
creating a pair of unpinned pancake-antipancake vortices fTOTd) . The new pancake follows the 

4 If the YBCO film is grown on a normal substrate, the transport current flows parallel to the ab planes 
and there is no breaking phenomenon. For this reason, it is required a miscut angle in the substrate in order 
to obtain a certain component in the c direction and observe vortex breaking. 
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whole vortex drift, while the antivortex approaches to the original pinned pancake, anihilating 
to each other. Finally, we obtain a depinned vortex with undeformed SV segments, Fig. [TUlfe). 
for a driving force much lower than that one required for pancake depinning. 
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Figure 10: Static vortex deformation (a,b) and proposed breaking and cross-joining process (c,d,e) 
for a string vortex (SV) in a vicinal film with a field orientation 8 = 10°, Fig. [U The graph is the 
projection in the ab plane of a kinked vortex, where lines are SV between a-b planes and circles are 
pancake vortices. The solid (red) lines are for SVs in the graph plane, dash (maroon) ones are for 
SV at lower planes and dot (orange) lines are SV in higher planes. Solid and hollow dots are pancake 
vortices and antivortices, respectively. 



This vortex breaking and cross-joining process is set up when the driving line force is larger 
than the sum of the maximum elastic and pinning line forces. Using the angular coordinates 
of Fig. (9[ the driving line force in a SV is = &oJ sin^fo, where 9 V is the vicinal angle and 
b is the unit vector in the b axis. From Eq. (]18p . we obtain that the maximum elastic line 
force is f e = — (2e^ / c?ii) | tan#|b, where d\\ is the a-b interlayer distance. Considering this and 
a certain pinning force, we obtain the following critical current density 

J c = |tanfl| + -^-, (45) 

consistent with the expression in Refs. [22j [28] with a total SV breaking force of 2e;. 

5.2.2 Comparison with experiments 

The model introduced above can be applied in order to obtain e; for a SV, as follows. In Fig. 
[TT]we present the measured J c (9) for an YBCO vicinal film with 4° miscut substrate [28J. In 
this figure we include the fitting curve for the breaking limited J c (9) of Eq. (|45p (solid line) 
together with the fitted J c (9) for a kinked vortex depinnig (dash line), calculated using [28] 



T _ jp.str cos 9 + f P:PC sin 1 9 \ 

Jc ~ $ 8m\e - e v \ ' 1 j 

where / PjS tr and /p, pc are the pinning force per unit length for the string and pancake vortices, 
respectively, and 9 V is the substrate miscut angle. The measured data of Fig. [11] is optimally 
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Figure 11: Measured J c {9) for an YBCO vicinal film with 4° miscut substrate [28] (line with circles) 
together with the fitting curve for the breaking limited J c of Eq. (|43|) (solid line) and the fit for a kinked 
vortex depinnig (dash line), calculated using Eq. (j4"6")l . The fitting parameters are ei = 1.69 x 10~ 14 N, 
f p<atr = 6.62 x 10- 6 N/m and f p . pc = 1.82 x 10~ 4 N/m. 



fitted by Eqs. (@5J) and (HHD with e, = 1.69 x 1(T 14 N, / PjStr = 6.62 x l(r 6 N/m and / PiPC = 
1.82 x 10 _4 N/m. The two vertical dot lines in Fig. [TT] delimit the expected angular region 
where kinked vortices can exist {0 ~ 11°) [3J, while the dash-dot line indicates the 9 value at 
which the kinked vortex presents null total driving force, 9 = V , so that the depinning- limited 
J c diverges. Thus, it is clear that for this situation the limiting mechanism of J c must be 
flux breaking. Besides, both Eqs. (|45l) and (|46p for string vortex breaking and kinked vortex 
depinning, respectively, fit to the experimental data for the range \9\ < 9°; a slightly smaller 
value than the originally expected \9\ ~ 11°. 

The expected SV line tension in the parallel direction from Eqs. (jlip and (|44p is e\ ~ 
1.9 x 10 _12 N, which is around two orders of magnitude larger than the obtained one from 
Fig. HU This discrepancy, also found in Sec. 15.11 evidences that the use of the simplified 
Eqs. ©-((H]) is probably too naive and a more detailed study of the line tension in terms of 
the flux line lattice elasticity is required for a fine prediction. 

6 Conclusions 

One way of studying the superconducting vortex lattice properties is by vortex breaking 
experiments in nonuniform pinning force distributions. 

In this article, the vortex bending and breaking in a high-pinning/low-pinning/high- 
pinning layered structure has been studied from a microscopical point of view in order to 
extract vortex internal properties from breaking experiments. Our study has been based on 
an elastic model assuming a constant vortex line tension, ej. For the general 3D bending, this 
model is applicable for isotropic superconductors in the single-vortex approximation. The 
latter approximation can be done for a magnetic penetration depth A much larger than the 
vortex separation because for this case the vortex line energy is mainly the core energy. As 
a consequence, the line tension is independent on both the vortex deformation amplitude 
and the deformation length. Within the assumption of constant line tension, the model is 
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applicable to any nonuniform pinning force or current density distribution. 

Summarizing, the main results from the elastic model are the following. 

When the pinning force in the high-pinning regions (HPRs) is much higher than that in the 
low-pinning ones (LPR), the vortex only bends in the LPR. For low magnetic field orientation 
angles with the HPR-LPR boundary (ip), the vortex takes the shape of a circular arch with a 
minimum curvature radius of half the LPR thickness. The latter situation corresponds to the 
maximum elastic force. For a transport current J higher than the value which this happens 
(J c ), the elastic force cannot balance the driving one and the vortex moves until it finds its 
neighbour. After this, we have proposed that the adjacent vortices cross-join and recombine 
with an steady average vortex movement in the LPR. For higher B orientation angles \<p\, 
there is a similar behaviour, with the difference that the vortex deforms in two ortogonal 
directions. When <p is close to 7r/2, corresponding to the force- free configuration, the vortex 
experiences a helical deformation with a small radius. For the case that ip is exactly vr/2, the 
vortex presents a helical instability for J above a certain threshold and keeps as a straight 
line for J below. 

The more realistic case of finite pinning force in the HPR presents vortex bending in both 
in the LPR and the HPRs. For this case, if the maximum elastic line force in the LPR is 
larger than the pinning line force difference, the vortices present depinning in the HPR above 
a certain orientation angle ip^. When the orientation angle is slightly lower than the 
vortex bends in the HPR over a thickness larger than that of the LPR. 

The model calculations have been used for obtaining the vortex line tension in two ex- 
perimental situations: a low-angle grain boundary (LAGB) in an YBCO bicrystal and an 
YBCO film on a miscut substrate. In these systems the vortices involved are vortices of 
Abrikosov-Josephson nature and string vortices. 

For LAGB with a low angle <p between the magnetic field and the grain boundary plane, 
the vortex mainly bends in one single direction and the constant line tension assumption is 
reasonable. The line tension obtained by fitting the modelled curve with the measured data 
is much smaller than the expected one for an anisotropic vortex from Eq. (jlOp . although e; is 
of the same order of magnitude as an isotropic vortex with penetration depth A c . This result 
is consistent with an anisotropy loss due to a vortex enlargement in the c crystallographic 
axis. However, it must be kept in mind that the applicability of the simplified equations 
(]10p - (llip is not evident for the studied LAGB situation. 

The results for YBCO films with a miscut angle show again that the line tension is two 
orders of magnitude lower than expected. This reinforces the hypothesis that Eqs. (I10p - (|lip 
are not applicable for SVs at the studied temperature and magnetic field conditions, for which 
the vortices are strongly overlapped. Thus, a more complex analysis in terms of the flux line 
lattice elasticity should be done for performed for a proper prediction of the line tension. 

In conclusion, modelling based on the vortex elastic properties provides a good explanation 
for the flux breaking phenomena. In this article, the elastic theory has been successfully 
applied for obtaining the line tension of "exotic" vortices, such as string vortices between ab 
planes and in low-angle grain boundaries. 
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